jordan_sqrt <-
function(x){
	out <- sqrt(lambda1(x))*u1(x) + sqrt(lambda2(x))*u2(x)
	return(out)
}

